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Abstract 

Radioactive decays contribute significantly to the re-heating of supernova ejecta. Previous 
works mainly considered the energy deposited by 7-rays and positrons produced by ^^Ni, ^^Co, 
^'^Ni, ^'^Co, ^^Ti, and ^^Sc. We point out that Auger and internal conversion electrons constitute 
an additional heat source. At late times, these electrons can contribute significantly to supernova 
light curves for reasonable nucleosynthetic yields. In particular, the internal conversion electrons 
emitted in the decay of ^^Co are an important heating channel for supernovae that have become 
largely transparent to 7-rays. We show that when the heating by these electrons is accounted for, 
the slowing down of the light curves of SN 1998bw and SN 2003hv is naturally obtained for typical 
nucleosynthetic yields. Additionally, we show that for SN 1987A the effects of internal conversion 
electrons are likely significant for the derivation of ^'^Ti yields from its late time bolometric light 
curve. 



1 Introduction 

A great success of nuclear astrophysics was the demonstration that radioactivity powers the hght curves 
of Type la supernovae (SN la, most hkely thermonuclear) and, at least for some events at late times. 
Type Ib/c and Type II (most likely core collapse) supernovae. Shortly after its explosion a supernova 
enters a state of homologous expansion. The temperature drops during the expansion until nuclear 
fusion reactions cease; radioactive decays, however, still take place, since the matter is only slightly 
ionized. It is now widely accepted that the energy liberated in the decay chain of radioactive ^^Ni is 
the most important nuclear source for re-heating the supernova ejecta to temperatures high enough for 
the spectrum to peak at optical wavelengths [H |2]. At first, the bulk of the heating is produced by the 
energetic 7-rays which thermalize and deposit their energy via Compton scattering and photoelectric 
absorption. In the homologous expansion, the column density (and therefore also approximately the 
Compton opacity) decreases with time as and the ejecta become more and more transparent to 
these high energy photons. Once 7-rays escape, the positrons produced in the decay of ^^Co and '''^Sc 
were thought to be the main heating sources. Here, we draw attention to often overlooked additional 
leptonic heating channels: Auger and internal conversion electrons. In section [2] we review the physics 
of nuclear decays relevant for supernova light curves. In section [3] we demonstrate the impact of internal 
conversion electrons from the decay of ^''Co on different supernova light curves. We conclude in section 
Hlwith an outlook how this effect could be used to constrain supernova explosion models in the future. 
For the published refereed journal article that first pointed out the significance of internal conversion 
electrons on supernova light curves, please see [3]. 
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2 Nuclear decays 



The time dependence of n nuclide abundances Ni in a ecay chain is governed by the Bateman equations: 

^ = -''"^ 

^ = \.,N,.,-\,Ni . (2) 

at 

For n = 2 and initial abundances A''i(0) and A''2(0) we get the solution 

N,{t) = iVi(0)exp(-Ait) (3) 

N2it) = iVi(0)-^^[exp(-Ait) -exp(-A2t)] + Ar2(0)exp(-A2t) . (4) 
A2 — Ai 

1/2 

The decay constants Aj are related to the half-lives and the mean life-time Tj via 

The rate of energy deposition by decays of nucleus i is given by the activity multiplied by the energy 
deposited per decay: 

ei = \N,{t)qi{t) , (6) 

where the number A'^j is given by eq. Oor eq. [Hand the energy deposited, Qi, is a function of time due to the 
increasing escape fraction of 7-rays and possible late time escape of positrons. 

Most of the nuclei synthesized in supernova explo- Table 1: Radioactive decay energies (keV decay"^) 
sions have atomic masses in the range A ^ 12 — 70. 
The isotopes are generally either stable or on the 
proton-rich side of the valley of stability. For these 
unstable nuclei, radioactive decay occurs along an 
isobar towards neutron richness, either by electron 
capture or positron emission. Electron capture pro- 
ceeds via the capture of an atomic (typically K or L 
shell) electron by a proton in the nucleus and corresponding emission of an electron neutrino. Positron 
emission proceeds via the decay of a proton in the nucleus into a neutron and the corresponding emis- 
sion of a positron and an electron neutrino. For both processes, the transition to the daughter nucleus 
is a statistical process to a distribution of (excited) nuclear states. Following electron capture, the 
daughter is formed with a hole in its atomic electron cloud. Higher lying atomic electrons transition to 
fill the gaps in the lower lying atomic shells, which results in characteristic X-rays being emitted from 
the cascade. For every such electron transition, there is also a probability that instead of an X-ray 
one or more higher lying atomic electrons are ejected. These electrons are known as Auger electrons. 
Analogous to the X-ray cascade of the atom, the excited daughter nucleus typically undergoes a series 
of transitions towards the ground state accompanied by emission of characteristic 7-ray photons. If the 
nucleus is surrounded by an electron cloud, then there is again a probability that instead of a 7-ray one 
or more atomic electrons are ejected. This process, in which the energy difference of the nuclear levels 
is carried away by the ejection of an inner atomic electron and there is no 7-ray photon emitted, is 
called internal conversion. The probability for internal conversion to occur for a given level is measured 
by the internal conversion coefficient a = dc-cxcitations ^ ^ normally small for nuclei with low 

#01 7 de— excitations 

atomic number, but increases for levels close to the ground state. For the decay of ^^Co the probability 
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for production of internal conversion electrons is very large, due to the existence of a low-lying nuclear 
level in the daughter nucleus ^^Fe (see Fig. [T]). For the first exited state (14.4 keV 3/2~) of ^''Fe, the 
internal conversion coefficient a = 8.58. This decay is 100% electron capture, and over 99.8% all decays 
are into the 136 keV level. Captures to the higher lying 366.74 and 706.42 keV levels contribute only 
marginally. For ground-state to ground-state electron capture transitions, such as in the decay of ^^Fe, 
no 7-rays or positrons are emitted. Assuming that the neutrino escapes without interactions, in such 
transitions the Auger electrons and the X-rays constitute the only sources of radioactive heating. 
The following four decay chains contribute most to bolometric 
supernova light curves: 
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Here we do not model the radiative transport, we only compare 
leptonic (and X-ray) energy injection rates. We do not include 
the energy produced by the pair annihilation and further assume 
that the kinetic energy of the leptons is completely and locally 
deposited and thermalized. The energy generation rates presented 
here do not include any heating due to 7-rays and therefore are not 
predictions for bolometric light curves. This approach, however, 
does allow for a direct comparison of the relative importance of 
the different leptonic heating channels (positrons from the decays 
of ^^Co and ^^Sc and electrons from the decays of ^^Co and ^^Fe.) 
The relevant energies of the different decay channels are listed in 




table [H These data are taken from the National Nuclear Data Centei0. 



14.41 keV 



Figure 1: Ground state and first 
two exited levels of ^^Fe showing 7- 
ray energies and intensities for the 
decay of ^^Co. 



3 Late-time bolometric light curves 

Bolometric light curves have to be reconstructed from multi-band photometry. For reliable reconstruc- 
tion, the contribution of the UV/optical {UBVRI) and near-infrared {JHK) bands have to be included 
(UVOIR light curve). However, near-IR observations are rare at very late epochs, and sometimes only 
5-through-/ band observations with a near-IR correction extrapolated from earlier epochs are used [1] . 
Below, we discuss some aspects and examples of different supernova light curves. 

3.1 Type la supernovae 

SNe la are thought to be thermonuclear disruptions of white dwarf stars and as such do not have an 
extended envelope. Consequently, their ejecta become transparent to 7-rays relatively soon and the 
positron-dominated phase generally starts 150-300 days after the explosion |5l [6]. Between 300 and 
600 days, the UVOIR light curves fall exponentially with the ^^Co decay half-life [7], a clear sign that 
a constant fraction of positrons is trapped. Unfortunately, only few SN la were observed at even later 
times, but there are indications of a slow-down in the light curves of some optical bands after ~ 600 
days O [8]. A particularly interesting case is SN 2003hv, which shows a slow-down in the bolometric 

^http:/ /www. nndc.bnl.gov/ 
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Figure 2: Instantaneous energy generation 
rates for initial abundances taken from the 
W7 model of [11]. Thick lines are due to 
positrons, IC and Auger electrons alone. The 
thin lines also include the full X-ray dose, 
which is small even in the limit of complete 
X-ray trapping. The contribution of elec- 
trons and positrons from the decays of ^'^Ti 
and ^^Sc is too small to be seen in this figure. 
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light curve 786 days after B-band maximum light [9]. It has been suggested that additional (possibly 
radioactive) heating sources may be needed to explain this observed effect. We point out that the 
magnitude and time of occurrence of the observed slow down of the light curve is expected and a 
natural consequence of the leptonic energy injection from the decay of ^''Co. To illustrate this, consider 
the leptonic energy generation rates of important long lived isotopes between 500 and 2000 days for W7 
[To], a common reference SN la model (see Fig. [2]). For this particular choice of yields |TT], the light 
curve is significantly higher from assuming only heating from ^®Co positrons after about 750 days. 1000 
days after the explosion, the heating from ^^Co electrons equals the heating from ^^Co positrons. 

3.2 Core collapse supernovae 

The progenitors of core collapse SNe are massive stars, which have an extended, massive envelope at 
the time of explosion, which provides opacity for the emerging 7-rays. Compared to SNe la, the nuclear 
burning also takes place at higher entropy, which means that core collapse supernovae synthesize more 
44rp. ^j^ggg (differences imply that the moment when internal conversion electrons from ^''Co significantly 
contribute to the light curves is delayed. Nevertheless, two observed examples of core collapse SN where 
internal conversion electrons do significantly alter the shape of the light curve are given below. 

3.2.1 SN 1998bw 



Figure 3: Instantaneous energy generation 
rates due to radioactivity for initial abun- 
dances taken from a hypernova model of [12] . 
Only leptons and X-rays (which are small 
compared to the leptonic contribution) are 
included (i.e. no 7-ray contribution or freeze- 
out effects are considered). Squares are the 
data of 1998bw of [1]. Arrows represent 3a 
upper limits. There are no free parameters. 
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The well observed SN IQQSbw, which was the first case of a 7-ray burst associated with a SN[T3], was 
classified as a SN Ic. The explosion was a very energetic , asymmetric disruption of a massive, stripped 
stellar core (> 10 Mq; [H]). The high ejecta velocities and the lack of an extended envelope reduce the 
7-ray opacity such that by 1000 days almost all 7-rays are freely streaming. A nucleosynthetic yield 
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calculation for hypernovae like 1998bw [15] predicts a ^^Ni to ^^Ni mass ratio n^V^^ ^ 0.0366, ~ 1.5 
times the solar value for the ^^Fe to ^^Fe ratio (which is 0.0234; [ISj). After modeling the UVOIR light 
curve to ~ 1000 days after the explosion, [1] showed that a simple model without contributions from 
freeze-out effects, circumstellar interaction, accretion onto a central compact object or light echoes 
requires a ^^Ni to ^^Ni ratio ~ 13.5 times greater than solar. This large discrepancy between the 
value of 7?.^'^/s6 predicted from explosive nucleosynthesis calculations and the one derived from light 
curve modeling reduces if the effects of internal conversion and Auger electrons are included in the 
light-curve calculations. In Fig. |3] we show the combined leptonic and X-ray luminosity corresponding 
to the nucleosynthetic yield calculations of jl5]. The non-7-ray heating due to ^^Co is the dominant 



contribution between ~ 1000 and 1600 days. The observed slow down of the light curve of 1998bw at 
~ 900 days (see Fig. [3]) is thereby naturally obtained without the need for strongly super solar 7^^^/^^ 



3.2.2 SN 1987A 
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Figure 4: Leptonic energy deposition rates 
for important isotopes of the predicted yields 
for 1987A of [2l]. Note the hitherto ignored 
leptonic contribution of ^^Co (dot dashed) 
and its impact on the leptonic heating rate 
between 1000 and 2000 days. 



SN 1987A has demonstrated that light curves can be constructed for nearby supernovae extending for 
several years after the explosion. This has led to the consideration of longer lived radioactive species. 
The most important of these is ^''Co, which is expected to be produced in significant amounts as the 
decay product of the short lived ^^Ni. Due to its relatively longer half-life (271.79 days) and the high 
opacity of the ejecta to the emitted 7-rays, its 7-rays alone may dominate the bolometric light curve, 
especially of core collapse supernovae, at late times [T71|18]. In fact, due to its extended envelope, the 
opacity to 7-rays in the remnant of 1987A remained high for such a long time that it never entered the 
^^Co positron dominated phase [19]. It is argued, that the only other radionuclide which noticeably 
contributes to the bolometric light curve is ^^Ti [2U]. ^^Ti has an even longer half-life of 59.8 years, and 
the short lived daughter ^^Sc has a strong positron channel, which means that the ^^Ti decay chain is 
the dominant energy source at very late times. Photometric data of SN 1987A lead to the conclusion 
that the 7^^^/^*^ was rather large, about five times the corresponding solar value [211 122]. Due to the 
proximity of SN 1987A, the presence of ^^Co and ^^Co was inferred not only photometrically from their 
imprint on the light curve, but also through a direct detection of escaping 7-rays from their decays. The 
observed gamma ray flux from ^^Co with the Oriented Scintillation Spectrometer Experiment (OSSE) 
on board the Compton Gamma Ray Observatory [23] favored at most moderate enhancement of ^''Ni 
(factor 1.5), in good agreement with nuclear reaction network calculations [21]. Similarly, based on the 
light curve, a rather high ^^Ti abundance was claimed in several studies [22 ] [25 ] 126] . The great difficulty 
of explaining the light curve of 1987A between 1000 and 1500 days without strong overproduction of 
^^Ni can be ameliorated by including "freeze-out effects" [27] or by conveniently re-arranging the spatial 
distribution of the radionuclides [28], but some disagreement of data and models persists. By comparing 
the various leptonic energy deposition rates for a nucleosynthesis model of 1987A [23], we show that the 
leptonic channels of ^^Co decay are significantly contributing between 1000 and 2000 days (see Fig. H]). 
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4 Conclusion and Outlook 



Fitting models to observed late time light curves provides a unique and independent method to directly 
measure the isotopic yields of prominent radioactive nuclei synthesized in the explosion (in particular 
^^Ni, ^^Ni, and ^^Ti). We have shown that at late times, when the ejecta have become largely transparent 
to 7-rays, the energy carried by Auger and internal conversion electrons may constitute a significant 
source of heating. These additional decay channels have to be considered for reliable isotopic abundance 
determinations from light curves. In particular, we have shown that a re-analysis of the bolometric light 
curve of 1987A (taking the hitherto unconsidered effect of internal conversion electrons into account) 
would likely yield significantly different (smaller) ^^Ti and ^'^Ni masses. The new derived values for ^''Ni 
and ^'^Ti would allow us to gain more insight into the explosion mechanism of core collapse supernovae. 
In particular, observationally driven inferences about the location of the mass cut can be made, which 
would give us valuable information about the physical processes separating neutron star and black 
hole formation from massive stars. Last but not least, SNe la are considered as the source of the 
positrons required to explain the Galactic 511 keV annihilation line observed by Integral/SPI [29]. The 
question whether enough positrons can escape the remnant remains unanswered. An understanding 
of the physics that underlies the light curves of SNe is a crucial step in order to constrain the escape 
fraction of positrons. 
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